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Abstract. We construct a finite element like scheme for fully non-linear integro- 
partial differential equations arising in optimal control of jump-processes. Spe- 
cial cases of these equations include optimal portfolio and option pricing equa- 
tions in Finance. The schemes are monotone and robust. We prove that they 
converge in very general situations, including degenerate equations, multiple 
dimensions, relatively low regularity of the data, and for most (if not all) types 
of jump-models used in Finance. In all cases we provide (probably optimal) 
error bounds. These bounds apply when grids are unstructured and integral 
terms are very singular, two features that are new or highly unusual in this 
setting. 



1. Introduction 

In this paper we introduce and analyze finite element (FEM) like schemes for 
nonlocal Hamilton- Jacobi-Bellman equations (HJB equations) of the form 

(1.1) 

m 



sup {— tr [a{x, v)D'^u] — b{x, v)Du + c{x, v)u — f{x, v) — X^u{x)} =0 ir 



or 



(1.2) 




sup { 


— tr [a{x, v] 






where 




(1.3) 


a(x, v) 


(1.4) 


ruix) 


(1.5) 


J^uix) 



[u{x + ri{x, V, z)) — u{x)]v{dz), 

[u{x + ri{x, V, z)) — u{x) — l|2|<ir7(x, v, z)Du{x)]h'{dz), 

E 

y is a compact metric space and E = M*^ \ {0}. The coefficients a,ri,b,c,f are 
loosely speaking Lipschitz continuous in x, continuous in v, and Borel measurable 
in z. The precise assumptions will be given in the next section. These equations 
are the dynamic programming equations for stochastic control problems involving 
Levy processes, a class of Markov processes with jumps [TT]. The two forms of the 
equation correspond to different intensities v{dz) of small jumps in the correspond- 
ing Levy processes. In general these equations are degenerate and fully-nonlinear, 
and the corresponding solutions are typically only Holder continuous and have to 
be understood in the viscosity sense [25l [291 H] ■ 
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Equations like and (|1.2p appear in advanced models of financial markets 

where price evolution of stocks (and other risky assets) are modeled as (exponential) 
pure-jump or jump-diffusion processes. Special cases of (jl.ip and (|1.2p include linear 
equations used in option pricing problems, obstacle problems used e.g. in pricing 
of American options, and fully non-linear equations (the full HJB equation) used in 
optimal portfolio problems. We refer to [19] for the first two cases and to [25j[8] for 
the last case. It is well-known that the standard Black-Scholes model (a diffusion 
model) give poor fit to real markets, at least on smaller time scales. E.g. log- 
returns distributions of stock prices are leptokurtic and have longer and fatter tails 
than predicted by the Black-Scholes model, see e.g. [12]. To improve upon these 
shortcomings, many pure-jump and jump-diffusion models have been suggested in 
the literature over the years, see e.g. [TOl [iP] for the most popular models. The 
empirical fact that Levy processes with discontinuous sample paths tend to better 
model e.g. stock prices, is one main reason for the renewed interest in stochastic 
optimal control of jump-diffusion processes. 

Except in very simple cases equations (jl.ip and (|1.2p do not have closed form 
solutions, so numerical methods are needed to obtain (approximate) solutions. In 
this paper we construct FEM like schemes. Since the equations are fully nonlinear 
we have no weak/ variational formulation, and the usual FEM approach does not 
work. If the measure v is finite and the corresponding HJB equation is (|l.ip . 
we discretize in two steps. First we obtain a semi-discrete equation ("discrete in 
time") as the dynamic programming equation of a discrete time control problem 
approximating the underlying continuous time control problems of (|l.ip : 

(1.6) 

Uhix) = inf <^ h](x, v) + e-'*'^'^''') — — V [uhix + hh{x, v) -f Vham{x, v)) + 

m— 1 

+ Uh{x ^ hh{x,v) — \fhara{'X,v)^ ^ ^ y Uhix ^ r\(x,v,z))v{dz) |, 

where Om is the m-th column of the matrix cr, ft, is the discretization step and A is the 
mass of the measure. The next step is obtain a fully discrete equation by introducing 
a regular triangulation and look for continuous piecewise linear functions over the 
chosen triangulation satisfying the semi-discrete equation at every vertex of the 
triangulation. When v is not finite we first approximate it by a finite measure 
^T<\z\<By{dz) (truncation), and then approximate the truncated equation following 
the above approach. To improve the truncation approximation, we also add small 
diffusion and/or drift terms to the equation. 

We prove that these methods converge and derive (probably optimal) error 
bounds using the framework developed in [2D] and ideas from [T7] . We also discuss 
issues like restricting to a bounded domain, truncating long jumps, and approxi- 
mating integrals by quadrature. In all cases we provide rigorous error bounds for 
the various approximations. What remains to do in a computer implementation, is 
the resolution of the non-linearity. There are various ways to do that, e.g. policy 
iteration, value iteration, artificial time methods (stationary limits), and so on. We 
will not address this point in this paper. 

FEM schemes like the one we have described above are usually called semi- 
Lagrangian or control schemes in the literature. They are (usually) monotone and 
first order accurate, and a comprehensive background and references can be found in 
[21] . Most results in this field concerns deterministic control problems and first order 
HJB equations without integral terms. Semi-Lagrangian schemes for second order 
HJB equations with no integral terms have been considered in j37l HTl [S] [1] . More- 
over in ^16j such schemes were derived for HJB equations associated to piecewise 
deterministic processes (compound Poisson processes with drift). One advantage 
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of this type of schemes is that in general they produce also an approximation of 
the optimal control law in feedback form and the optimal trajectories (this is a key 
point in the study of a control problem). This advantage can also be shared by some 
(monotone, low order) finite difference schemes following the construction of Kush- 
ner, and we refer to [34] for a discussion of this point. An other advantage is that 
these schemes are formulated on general grids/triangulations just like FEMs. By 
contrast this is very cumbersome to achieve with finite difference methods. A third 
advantage is that the semi-Lagrangian scheme automatically handles non-diagonally 
dominant diffusion matrices a. Such matrices appear in applications and are cum- 
bersome to discretize, we refer to [12] for a finite difference (FDM) approach to this 
problem. 

The construction and analysis of numerical schemes for linear integro-partial dif- 
ferential equations arising as pricing equations in financial markets of jump-diffusion 
type is currently an active field of research, see e.g. [HI [36l |23l [15] and references 
therein. By contrast, there are few works on numerical schemes for fully nonlinear 
degenerate integro-partial differential equations. We mention the discussion about 
jump-diffusion processes in ^34] and the papers [221 [HI [30l [10] . In cases mainly 
monotone FDMs are considered and convergence is obtained. In the last two papers 
the main focus is on convergence rates. The framework of [30] is again based on 
ideas of Krylov, Barles and Jakobsen [311 [32l [5] in the pure PDE case. We em- 
phasize that the error bounds of this paper apply when grids are unstructured and 
integral terms are very singular, a feature which is new or highly unusual in this 
setting. 

The rest of this paper is organized as follows. In section [5] we state the assump- 
tions on the data and give well-posedness/regularity results for equations (jl.ip and 
(|1.2p . We discuss truncation of the Levy measure and reduction to a bounded do- 
main and give error bounds. In section [3] we construct the schemes via the dynamic 
programming equation (a semi-discretization) and FEM ideas. We prove existence, 
uniqueness, consistency, and partial error bounds. We also derive a fully discrete 
scheme of FEM type by piecewise linear reconstruction on a (regular) triangulation 
of the domain. In section|4]we derive error bounds for the semi-discrete scheme with 
or without truncation of the Levy measure. Finally in the appendix, we prove the 
main technical result of this paper, a regularity and continuous dependence result 
for the semi-discrete scheme. 

Notation: By [/S'Cb(R^) and Cfc(R^) we mean the spaces of bounded continuous 
and upper semicontinuous functions. We will use the following norms 

l/lo = esssup^gRiv|/(2:)|, [f]i = \Df\o, and |/|i = |/|o + [./li- 

2. Preliminaries 

In this section we state the assumptions for equations (|l.ip and (II. 2p and give 
well-posedness results. We discuss reduction to a bounded domain and reduction 
to bounded measure v with compact support. 

2.1. Assumptions, well-posedness, regularity. We will use the following as- 
sumptions: 

(Al) The set y is a compact metric space, the coefficients a,r],b,c,f are con- 
tinuous in X and v and Borel measurable in z, and ^ is a positive Radon 
measure on E := M^^ \ {0}. 
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(A2) There exists Li, L2,l > such that for any v ^ V and z E E 

\ci;v)\i + \fi;v)\i<L,, 

\a{-,v)U + \b{-,v)\i<L2 

<i2(k|l|.|<i+e^l^ll|,|>i). 
(A3) There exists cq > such that, for all x E and v E V, 

c{x, v) > Cq. 

(A4) There exists L3 > 0, such that for £ defined in (A2), any v eV and x eR^, 

\DMx,v,-)\o<L3e'\'^. 

For the measure v we will use the following integrability assumptions. 
(Bl)fe The measure v satisfies for £ defined in (A2) and some k E {0, 1, 2}, 

/ (|z|'=l|,|<i+e^l^ll|,|>i)^('^^) <oo. 
Je 

(B2) The measure i' has a positive density m : E [0, 00) such that i^{dz) = 
ra{z)dz and for some C, £ > 0, and n e N 

\D''m{z)\<Cke-^^+'^\'\ for fc = 0,l,...,n, 

where D^m is the vector of all order k derivatives of m and £ is defined in 
(A2). 

(B3) The measure v has a positive density m : ^ [0, 00) such that v{dz) = 
m{z)dz and for some a E [0, 2), Co, Ci, and e > 0, 

\D^m{z)\ < Ck{ ^^^^\^^^ l|.|<i + e-(^+-)l-ll|,|>i) for A: = 0, 1, 

where D'^m = m, D^m — Dm the gradient of m, and £ is defined in (A2). 
Assumptions (Al) - (A3) and (Bl)2 are standard from a stochastic control theory 
point of view, as they insure existence and uniqueness of strong solutions of the 
underlying stochastic differential equations, see [26]. Under assumption (Bl)2 or the 
less general assumption (Bl)i the measure v may have a (non-integrable) singularity 
at z = 0. If (Bl)i or (Bl)o holds the HJB equation takes the form ((LT]). If (Bl)2 
holds but not (Bl)i, then (|1.2p gives the correct form of the HJB equation. In 
this case, the extra term in the integrand of J^^ (compared to T"") is needed for 
the integral to converge: If (Bl)2 holds, then J'"u converge for all functions u 
is with polynomial growth at infinity. Finally, (B2) and (B3) prescribe densities. 
Close to z = 0, the density of (B3) equals the density of the a-stable processes 
related to the fractional Laplacian A"/^. (B3) implies (Bl)2, and if a S [0, 1), then 
(B3) also imphes (Bl)i. 

Example 2.1. When the stock returns are modeled as exponential Levy pro- 
cesses the integral terms equals (jl.4p or ()1.5|) with ri{x, z) = x{e^ — 1) in the 
one-dimensional uncontrolled case, see [19;. An important extension of this ex- 
pression satisfying (A2) and (A4) is 77(0;, v, z) = fj{x, v)(j){z) where f] is bounded and 
Lipschitz and \(t>{z)\ < C(el^l - 1), \D(l3{z)\ < Cel^L 

Example 2.2. Assumptions (B2) and (B3) is satisfied by all bounded and un- 
bounded Levy processes used in the literature to model financial markets, see [TO] 
for a nice overview. In the models of Merton and Kou, the Levy measures have 
bounded densities (no singularity), in the Merton case [35] given by 

A I 1^ 

v{dz) — — dz for constants X^S,fi. 
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For these models (B2) holds. The Variance Gamma model has a Levy density with 
an integrable singularity at 2; = 0, and satisfy (B3) with a = 0. The Normal Inverse 
Gaussian model ^40^ has a non-integrable density corresponding to a = 1 in (B3), 
and the Levy measure is 

C 

v{dz) = —e^'' Ki{B\z\) for constants A, B, C, 
\^\ 

where Ki is a modified Bessel function of 2nd kind. Finally we mention models 
using tempered a-stable processes, e.g. the CGMY model Here a G (0,2) in 
(B3) and the Levy measure is give by 



Lz>0 



for constants a e (0, 2), c_, c+, A_, A+ > 0. 

It is well known that under the above assumptions the solutions to (|l.ip and 
(|1.2|) need not be smooth, and that the correct concept of (weak) solutions is that 
of viscosity solutions. For the definition of viscosity solution in this case we refer 
to [U 131 [Ml m ■ We state without proof a well-posedness and regularity result 
for (|1.2p . The proof of this result is standard, and we refer to [31 [351 UHl US] for the 
proofs of similar results. 

Theorem 2.1. Assume (Al) - (A3) and (81)2 hold. 

(i) There exists a unique viscosity solution u G Cb(M^) of equation p.2p which 
is Holder continuous, i.e., there is a S € (0, 1] such that 



\u{x) — u{y)\ < C\x — for all x,y E 



(ii) There exists a constant ci > depending only on sup„[(T(-, i;)]i, sup^,[5(-, 
and 



s\ypy j^[rj{- ,v, z)Wv{dz) such that if c^ > Ci, then the viscosity solution u of ^ 



is Lipschitz continuous (6=1 above). 

(Hi) Let u,—v G [/S'Ch(M^). If u and v are respectively viscosity sub- and 
supersolutions of (jl.2p . then u < v in . 

Remark 2.2. Note that this result also holds for (|l.ip under assumptions (Al) - 
(A3) and either one of (Bl)i or (Bl)o. Simply write this equation in the form (|1.2p 
with b replaced by 6 = 6 — r)i^{dz). 

2.2. Reduction to bounded measure with compact support. Consider Levy 
measures v which are not bounded nor have compact support. We assume that 
(B3) hold and consider two cases: (i) a G [0, 1) and (ii) a G [1, 2). In the first case 
we truncate the Levy and compensate by adding a drift and a diffusion term. Let 
us introduce the "two-scales" truncated Levy measure 

(2.1) VrM{dz) := l,.<|2|<fli/(dz), 

where r G (0, 1) and R> 1. Clearly, Vr^R{dz) is a bounded and compactly supported 
measure satisfying assumption (Bl)o. When a G (0, 1) assumption (Bl)i holds and 
the HJB equation is (jl.ip . This equation is approximated by 

(2.2) 

sn^{~iv[a{x,v)D'^u] — b{x,v)Du + c{x,v)u — f{x,v)—T'"u{x)]=Q in R^, 
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where 

(2.3) a{x,v) = a{x,v) + — f ri{x,v, z)ri{x,v, z)'^ iy{dz) 

2 J0<|2|<r 

(2.4) b{x,v) — b{x,v) + / ri{x,v,z)v{dz), 

Jo<\z\<r 



(2.5) I^u{x) — / [(l){x + ri[x,v, z)) — (l)(x)] Vr.B.{dz). 

Je 

This gives a "third order" approximation of since 



|X> +{b- h)D4> + tr [{a - a)D'^4)] ~ X>| < iflD^- 



H-r,z)\'v{dz) 

Q<\z\<r 



In case (ii) we truncate the Levy measure and compensate by adding a diffusion 
term. When a G [1,2) assumption (Bl)i is not satisfied and the HJB equation is 
(|1.2p . We approximate this equation by 

(2.6) 

snp I— tT[d{x,v)D^u] — b{x,v)Du c{x,v)u — f{x,v) ~i^u{x)> — in M^, 
where d and is defined in (|2.3p and (|2.5p . and 

(2.7) b{x,v) = b{x,v) + l\^\^ir]{x,v,z)vr,Bidz). 

Je 

Again we obtain a "third order" approximation of the integral term {J^"" (j) this time) 
since 



|I> +{b- b)D(j> + tr [(a - a)D^<P] ~ J''(j3\ < K\D^<P\q 



H;;z)\',.{dz) 

Q<\z\<r 



Equations (|2.2p and (|2.6p are almost, but not quite, of the same form as p.ip . 
In [3D] it is proved that Theorem [O still holds for solutions of ([^ and The 
following result gives error bounds for the approximations (j2.2p . p.6p of equations 

(ini), (na). 

Lemma 2.3. ^sswme (Al) - (A3) and (B3) hold and let u,u,v,v € C°'^(R^) 
solve (II. ip . (12. 2p . (jl.2p . (j2.6p respectively. Then for r > small enough and R 
large enough we have 

|u-S|o + |«-z;|o <Ci(ri-"/3 + e-^^), 

/or some constant Ci independent of r and R. 

The result about v,v is proved in [3D], while the result about w, u is new but the 
proof is almost identical to the proof of the w, v result. We omit it. 

Remark 2.4. The compensating drift and diffusion terms are added to improve the 
convergence of the approximation (of the small jumps). The error is of order j-i-"/^ 
without the the compensating diff'usion (see [30j), while it is of order r^^" without 
compensating drift and diffusion in the case a G [0, 1). 

Remark 2.5. In probabilistic terms, the explanation for the improved convergence 
is that the small jumps of a Levy process can be approximated by a Brownian 
motion, we refer to e.g. [2 [19] for details. Moreover using probabilistic methods 
(Berry-Essen type estimates), it is possible to prove convergence of order r in some 
cases [5] . We refer to [ID] where such estimates are made explicit for linear parabolic 
equations with constant coefficients. 

Remark 2.6. When a = no improvement is obtained by adding drift and diffusion 
terms to the equation. If v has a bounded density, no truncation is needed close to 
z = 0. 
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2.3. Reduction to a bounded domain. Reduction to a bounded domain is key 
step in order to implement a numerical method. Following ideas of [17j . we restrict 
equations (|l.ip and (|1.2p to bounded domains by truncating the coefficients outside 
some large ball, 

Bj_={xe M" : |a;| < -} for ^ > 0. 

Let e C^(M^) be a cut-off function satisfying < < 1 and £,fj.{x) — 1 for 
X £ Bi_ and define 

(^1,= S.pi{x)(j{x,a), bf, = ^f,(x)b{x,a), ■q^{x,v, z) = £,y,{x)-q{x,v, z). 

Since 5^, cr^ and ry^ satisfies the same assumptions of a, h and rj for any /i, there 

exists a unique viscosity solution of the equation 

(2.8) 

sn-p {-iv[ay,{x,v)D'^u]~ b^{x,v)Du + c{x,v)u - f{x,v) - J1l{x)] = Q in 

where a^(a;, = i(T^(x, v)(tJ^{x, v) and J7^(-) is defined as in p.Sp with 77^ in place 
of 77. Since the coefficients 6^, and 77^ are zero outside of supp(^^), the solution 
Ufj_ is given by 

uJx) = mm — tor any x f. supp(4m). 

v^v c\x,v) 

The next result give a crude bound on the error that this cut-off procedure 
introduces. 

Lemma 2.7. Assume (Al), (A2), (81)2 hold, that (fr2|l and ([2^ satisfy the 
dynamic programming principle, and that u, solve (|f .2p . ()2.8p respectively. Then 
there exists a constant C such that 

\u{x) - uJx)\q <Cn^{l + \x\^) in Bi. 

We skip the proof since it is similar to the proof of Proposition 3.f in [17] when 
we are equipped with the moment estimates (3.2) and (3.4) of |39]. The dynamic 
programming principle has recently been extended to the current setting in [11] . 

Remark 2.8. As in Remark 12.21 we immediately get an analogous result for (jl.ip 
under assumption (Bl)i. 

Remark 2.9. If the equation is uniformly elliptic or has a non-degenerate singular 
integral term, then you expect estimates decaying exponentially as fi 00. We 
refer to [201 136] for such result in a linear one-dimensional setting. 

3. Construction of the scheme 

In subsections 13.11 and 13.21 we always assume that (Al) - (A3) and (Bl)o hold: 
The Levy measure 1/ is bounded and the HJB equation has the form (II. ip . We can 
always reduce to this case by a truncation. In subsection 13 .31 we discuss the general 
case. 

3.1. Semi-discretization. We introduce a control problem for which (|l.ip is the 
corresponding dynamic programming or HJB equation. We start by defining the 
controlled dynamics. Since the measure v is bounded, we can normalize and 
obtain a probability measure /i as follows 

IJ.{dz) — —v{dz) 
A 

where 

(3.1) \ = [ v{dz) 

J E 
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Now we consider a Markov process X^' evolving according the SDE 

(3.2) dXt^b{Xt^,Vt-)dt + aiXt^,Vt^)dWt+ [ v{Xt- ,Vt- , z)fL{dz,dt), 

J\z\>0 

where /I is a Poisson measure corresponding to a compound Poisson process with 
jump intensity A and jump distribution fi. The control v. belongs to V, the set 
of all progressively measurable processes with values in V . Since we assume that 
ly and hence /i are bounded, on any finite time interval Xt will only jump finitely 
many times with probability one. Between two jump times Ti and Ti+i the process 
diffuses according to the SDE 

(3.3) dXt = KXt- , vt- )dt + a{X^- , v^- )dWt. 

For vt = V <E V and a{x,v) — ^a{x,v)(j'^ {x^v), the infinitesimal generator of the 
process Xt is 

L-"^{x) = lini ^ tr [a{x, v)D^ib] + h{x, v)Di; + I^^P{x) 

for ip e C^(M^). On the paths of the process Xt we define the discounted cost 
functional 

poo 

Jix,v.)^E, / f{Xt,Vt)e-fo<''^^^''>'''dt 
'-Jo 

and we consider the corresponding value function 

u{x) = inf J{x, v.). 
V. ev 

In [11] it is proved that u is the unique viscosity solution of equation (|l.ip . 

Following the approach of [16| I17j we construct an approximation scheme for the 
equation (|1.1|) by discretizing the associated control problem. We fix a discretization 
step h > and consider two stochastic processes iV„ and Z„, n g N, taking values 
in N and in and representing the n-th jump time and the corresponding z-jump 
(size and direction) of the Poisson measure p.. We set A'o = and Zq — and 
assume that Nn has independent /lA-exponentially distributed increments, i.e. the 
probability distribution of iV„ is given by 

¥[Nn+l -Nn> j|iVo, iVi, . . . , Nn] = 6^''^^ n = 0, 1, 2, . . . , 

while the n € N, are i.i.d. random variables with probability density fi. Now 
we define a discrete time stochastic process Xn approximating the continuous time 
process Xt- 

' Xq ^ x, 

Xn = Xn-1 + hb{Xn^l,Vn-l) + VhJ2m=l O'm (^n- 1 , 1 )C™-1 1 

for n = A^, + 1, TV, + 2, . . . , iV,+i - 

for i = 1,2,3,..., where Cm denote the m-th column of a and m = 1, . . . , d are 
random variables taking values in {—1,0,1} such that 

P[{e = ±l}] = ^ and P[U; 0} n {e ^ 0}] = 0, z^j. 

The discrete control {vn} is a random variable with values in V which is measurable 
with respect to the ct- algebra generated by Xi, . . . , Xn- 

Between jumps the process evolves like a random walk approximating the SDE 
p.3p . and when the process jumps there is no diffusion/random walk. The generator 
of the discrete process is 
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for ?/> e C"(R^) and where 
1 

/CJJV' = ^ + /i&(a;, w) + ^/ha„l{x, v)) + + hh{x, v) — A/7icr„j(x, t;)) — 2'4){x) 

rn — l 

with X'" as in (|1.4p . Observe that at this level the space variable is not discretized, 
therefore the discrete process has the same jump distribution as the continuous 
process. 

On the paths of the discrete process we define the cost functional 

oc 

(3.5) J^(x,K}) =E,[^/ie-''5:rr„^^(x„.,)/(x„,«„)], 

71=0 

(with the convention X]i=o ~ corresponding value function for the 

discrete control problem 

(3.6) Uh{x) = ini J,i{x,{vn}). 

Now it is easy to see, at least formally ^ , that the following dynamic programming 
principle holds: 

p 

Uh{x) - inf [ V he-'' ^'^o «„) + e"'^^" o <^^'-^)u,,{Xp+i)\ , 

for any P £ N. Taking P = in the above equation and noting that Xq = x, we 
get 

Uh{x) = inf E4hfiXo,v) + e-''<''»'-^UhiX,)] 



(3.7) = inf \hf{x, v) + e-''=(-''') [\— V {u^ix + hb{x, v) + Vh<J^{x, v j) 

m— 1 

1 - e"'''*' f 1 -I 

+ Uh{x + hb{x,v) - Vham{x,v)) ^ / Uh{x + r]{x,v, z))i'{dz) |. 



m— 1 

A 

Rearranging the terms in the previous equation and dividing by ft.e"'"^^^'") we get 

(3.8) sup <^ -Lluhix) + "-—Uh{x) - e'''=(^^''V(x, i;) I = in R^, 

where L'^{-) is as in p.4p . We will talk about sub- and supersolutions of this 
equation, meaning that p.8p holds as an inequality with < and > respectively. For 
the scheme p.Sp we have the following easy properties: 

Proposition 3.1. Assume (Al) - (A3) and (B1)q. 

(i) If Uh and Vh are bounded sub- and supersolutions of (|3.8[) . then Uh < Vh in 

(ii) Any solution u/j of p.8p is bounded and satisfies for all h > 0, 

I I . SUp„ \fi-,v)\o , . ■ /AO) 

\Uh\o < , where cq is as m (A J). 

Co 

(Hi) There exists a unique bounded continuous function uu solving p.8|) . 

(iv) The solution of p.8p is (at least formally) equal to the value function p.6p . 

(v) ForO<h<land(f>e C"'(M^) satisfying \(j)\o H h |£''*0|o < oo, 

\L-^{x)-Lli^{x)\ < CM\D^^\o + \D^^\o + \D^cl>\o)+C2hX{{l+\ f H)l^'/'lo + |^'0|o), 
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where the constants C'l and C2 only depend on sup^, |(t|o,sup^ |6|o. 

Remark 3.2. For the truncation of an unbounded measure to converge as —> 0, 
we need A ^ 00 as /i ^ 0, while the scheme p.8|) converges only if both h ^ Q and 
A/i ^ by (v). The last condition means that the small jumps must be resolved in 
the grid. 



Proof. We work with the scheme in the equivalent form (|3.7p . Note that the scheme 
is monotone, it has positive coefficients. If u and v be sub- and supcrsolutions of 
(13. 8p . an easy and standard computation using (j3.7p and assumption (A3) shows 
that 

Uh{x) - Vh{x) < e-^°''\{uh - Vh)+\o, 



and hence (1 



-c„h 



)\{ufi — w/i)^|o < which proves (i). A similar computation 



shows (ii) after noting that ^hc_i < 1 for he > 0. Next denote the right hand side 
of (|3.7p by ThUh, and note that Th is contraction in the | • |o norm, 

ThUh - ThVh < e^'^''°\uh - Vh\o- 

Existence and uniqueness of a continuous bounded solution to (|3.8p follows from 
Banach's fixed point theorem and this proves (iii). Part (iv) follows from the from 
the dynamic programming principle as explained before the proposition. Finally, 
to prove (v), note that since = A~^I^ we may write 



L> - Lii, = (£- - cDij + {r-r)ij - (1 - 



-hX 



)ay^-il~hX-\l^e-'^^))r,p 



where C^tp^x) — tr [a(x,v)D'^ip{x)] + b{x,v)Dip{x). By Taylor expansion, e.g. 
\{C-£l)^\<h\b\l\D'yj\oHh'\b\l+h\b\o\a\l^^^^ 

and the estimates |1 — e^^|, |1 — a::^^(l — e^^)| < the result follows. □ 

Next we give an optimal Lipschitz regularity and continuous dependence on co- 
efficients result for the scheme 



Proposition 3.3. Let and he solutions of p.Sp corresponding to the data 
a,b,c, f^r],!/ and a,b,c, fji],!/ respectively, assume both sets of coefficients satisfy 
(Al) - (A3) and (B1)q, and that hX < Co and h e (0, 1]. Then there exist constants 
Ci,L,K > (only depending on the data and Co) such that if cq > c\ (cq as in 
(A3)), then for all h>0,x,ye M^, 



\uh{x) - Uh{y)\ < L\x -y\+ K sup 

vev 



|/-/|o + |c-£|o 



+ |6-&|o + |(t-ct|o+ / \r]{-,z)-fi{-,z)\^iy{dz) 



1/2- 



If the coefficients are equal and Uh = Uh, then this is a Lipschitz regularity result, 
while if X = y then this is a continuous dependence on the coefficients result. This 
result one of the main contributions of this paper, and it will play a key role in 
the next section where error bounds are derived. It extends similar results of [5] to 
equations with integral terms. The proof is rather long and technical and we have 
put the main bulk in the appendix. In the pure diffusion case, the current proof 
simplifies considerably the arguments of . 

Proof. The result follows from Theorem I A . II in the Appendix after writing 



rji' 



Ae 



Ae-^'» 

and noting that by the Cauchy-Schwartz inequality, 

1 - e-''^ 



-Xh 



riv, 



b-b- 



Ae 



-Xh 



< |6-6|o + (e'^^-l)A-i/2 / \,j-f,\^y 



1/2 




□ 
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Remark 3.4. When cq < Ci the solution to (|3.8p is only Holder continuous. We will 
not discuss this case here, but refer instead to 5 for how to obtain results in this 
case. 

3.2. The fully-discrete scheme. In this section we introduce a FEM like dis- 
cretization of (|3.8p yielding a fully discrete scheme. For a nice introduction to 
FEMs we refer to For fc > let T*^ — {S'j^}jgN be a non-degenerate triangula- 
tion of R^, i.e. a collection of iV-simplices S'^ such that 



U 5,* 



sup(diam S'j) < fc, pk < sup(diam i?gfc). 



where p <E (0, 1), diam denotes the diameter of the set, and Bgk is the greatest ball 
contained in Sj. We denote by X'^ = {xijigN the corresponding set of the vertices. 



and introduce the space of continuous piecewise linear functions on T'^, 

W'' ^{we C(R^) : Dw{x) is constant in S^}. 
Every element w in W'' can be expressed as 

w{x) = ^ f3i{x)w{xi), 

for basis functions (the so-called tent functions) (3i G W'^ satisfying (ii{xj) — Sij 
for i,j g N. It immediately follows that < Pi{x) < 1, X]ieN'^*('^) ~ 
compact support, and at any x G M.^ at most N + 1 j3iS are non-zero. The family 
{f3j}j is a partition of unity. On any simplex S^, the /?j's are called the barycentric 
coordinates of S^. 

The fully discrete scheme can then be formulated as follows: Find the function 
u £ that satisfies p.Sp at every vertex Xi £ ^ or equivalently, 
(3.9) 



L{xi) = inf |f 



j 



M,,{v)u{xj) + {l-e-^'')J2PvivHx-)\+hf{x„v)Y 

j 



for every Xi £ X^ . Here the matrices M{v) and P{v) are given by 



(3.10) 



M{v) 



1 



for -^ (u) = Pj{xi + hb{xi,v) ± ^/7i(T,„(x^, i;)) , and 



(3.11) 



Pij{v) = ^ y" (^ji^i +ri{xi,v,z))iy{dz). 



Note that M is a stochastic matrix, and for any in, only -I- 1 entries of any row 
of are non-zero. The matrix P is non-zero only if the vertex Xj belong to a 
simplex which has nonempty intersection with the set Xi + i]{xi, i;,supp(i/)) for all 
v£V. 

As a final step we also discretize Pij by (monotone) quadrature 



(3.12) 



Q 



Azm ■= 



{zj)ujj where 



Zj £ E, 



CJ, > 0, 



satisfying the error bound 



(3.13) EAz[(t>] 



<j){z)dz-QAM < K\\D(I3\\liAz for ^£W'^\K"). 



Here Az is the discretization parameter. All sensible tensor product quadratures 
satisfy this (first order) error bound, and the bound is optimal if the integrand (j) is 
not more regular. The monotonicity assumption luj > is satisfied for compound 
Gauss and Newton-Cotes types of quadratures in any space dimension and when 
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the order is less than 9 in the Newton-Cotes case. We refer to [30] for examples and 
a wider discussion of these issues. If we assuming that v has a density m satisfying 
(B2), we get the following final scheme: Find u € such that 
(3.14) 



dxi) — inf \ t 



-hc{xi,v) 



-\h 



where is defined in (|3.10p and 

Pij{v) = ^Qaz[I3o{x^ + ■q{xi,v,-))m{-)]. 

We have the following existence, uniqueness, and partial convergence result for 
(I3J41) . 

Theorem 3.5. Assume (Al) - (A4), (B2), ((XT^ . hold, hX < Cq, and 

h E (0,1]. Then there exists a unique bounded solution UhkAz € to (|3.14p . 
Furthermore, if the solution Uh of (j3.8p belongs to C^'^{MJ^), then 



\Uh - UhkAz\o < 



\uh\i 



2k + CAz 



Proof. Existence and uniqueness follow from a fixed point argument like in the 
proof of Proposition 13.11 To prove the error bound, note that since UhkAz (x) — 
Y.j&\Pj^^)uhkAz{xj), f3-i > 0, and J^Pj = 1: 

(3.15) 

\uh{x) - UhkAzix)\ <E (^f3j{x)\uh{x) - Uh{Xj)\ + Pj{x)\uh{Xj) - UhkAz{Xj)\ 

3 en 

< \uh\ik + Y I3j{x)\uh{xj) - UhkAz{xj)\. 

The last term can be estimated by using p.Sp and (|3.14p . Easy computations show 
that 



\uh{Xj) - UhkAz{Xj)\ 



< e 



— hcn ^ — Xh 



\uhkAz-Uh\a + e-'''»il-e-^'') 



-Xh\ 



'^PijluhiXj) - UhkAzix.j)\ 



^PijUhiXj) -^PijUhiXj) +^YP,jUh{Xj) - J J Uh{ 



x^ + r]{xi,z,v))u{dz] 



Note that ^ Pij = jQaz[1] = 1 by (|3.13p since J2 Pj — 1- Furthermore, 



A = 



< 



B = 



1 
A 

Az 



IkUh{xi + T]{xi,v, z))m{z)dz ~ Qaz IkUh{xi + 77(2;^, •))m(-) 

Az 



-^\\D[lkUh{xi+r^{xi,v,z))rn{z)]\\^^<^[\uh\iL:>.Co + \uh\oC-i^ j e 



{uh ~ IkUh){xi + ri{xi,v,z))v{dz) < k\uh\ 



ill) 



IE A 

where Ik(t>{x) — - (3i{x)(l){xi) is piecewise linear interpolation of (f> on X'^ and we 
have used (B2), (A4), and (|3.13p . Combining these estimates and (|3.15p . using the 
properties of Pi{x), and remembering that h e (0, 1] then gives the result. □ 



Remark 3.6. Since h € (0, 1], \ufi — UhkAz\o < 
the estimates obtained in [TB] in a different setting. 



C\uh\i k+Az 



h ' 



which is consistent with 
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Remark 3.7. If we use the cut-off procedure explained in Section r2.3l fa way of reduc- 
ing to a bounded domain), then equation (|3.14p gives a finite system of equations. 
In view of the integral term the new effective domain then becomes + supp(i/) 
(see Section [231) . 



3.3. Schemes for unbounded measures ly. We consider general unbounded 
Levy measures v under assumption (B3). There are two different cases: (i) a G [0, 1) 
with HJB equation (jl.ip . and (ii) a G [1,2) with HJB equation (|1.2|) . To derive our 
schemes, we first reduce to a bounded Levy measure Vr.R as explained in Section [2] 
The result are the approximate HJB equations (|2.2p and (|2.6p . These equations are 
then approximated by a slightly modified version of the semi-Lagrangian scheme 
p.Sp (or equivalently (|3.7p ) defined in Section [3T] 

We propose the following semi-Lagrangian scheme in case (i) 



(3.16) 



Vh{x) = inf \hf{x, v) + e^'^'^^^^^) V (vh{x + hb{x, v) + Vha+^,n{x, v)) 

vev I L 4a ^ — ' V 

m—l 

+ Vh{x + hb{x, v) + Vha^^m{x, v)) + Vh{x + hb{x, v) — Vha-+^m{x, v)) 

+ Vh{x + hb{x,v) - Vha-^rnix,v))) ^ / Vh{x + 'r]{x,v,z))iyr,R{dz) \, 

' K,R Je J 

and in case (ii) 
(3.17) 

Wh{x) = inf <^ hf{x, v) + e-'''=(^'") — — — V {whix + hb{x, v) + v^CT+,,„(a;, v)) 
vev L L 4a ^-^ 

m—l 

+ Wh{x + hb{x, v) + Vha-^mix, v)) + Wh{x + hb{x, v) — Vhd-+^jn{x, v)) 

+ Wh{x + hb{x,v) - Vha-^rnix,v)) ^ / Whix + ri{x,v, z))ur_R{dz) >, 

^r,R. Je J 

where Xr.R '■= i^r.^idz), b,b,Vr,R are defined in Section [21 and (T±.m is m-th 
column of 



(3.18) a±{x,v) = a{x,v) ± J / ri{x,v, z)ri{x,v, z)'^u{dz), 

y Jo<\z\<7- 

where the square root denotes the matrix square root. 

Remark 3.8. The additional terms in (|3.16p and ()3.17p compared with (|3.8p . enable 



us to use a i y J rjrf u instead of y cc"^ + / rjrf u as diffusion matrix. The consis- 
tency relation Proposition 13.11 (v) still holds, and if rj{x,v, z) = r]i(x,v)r]2{z), then 
the square root in (|3.18p equals Cr]i{x,v) where C is the precomputable constant 
matrix ^ J^^^^^^^r]2iz)T]^{z)i^{dz). 



Remark 3.9. These schemes are similar to p.8p . and can be derived in a similar 
way. The conclusions of Propositions 13.11 and 13.31 (when 6-terms are replaced by b- 
or 6-terms as defined in Section [2]) still hold for solutions of (|3.16p and (|3.17p . We 
refer to [3(T for the technical modifications needed to handle the integral term in 
the diffusion coefficients. 
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Remark 3.10. Previously bounded quantities may blow up as r — s- 0. Indeed by 
(B3) and (A2) we have for r e (0, 1), 

(3.19) 



K,R<^ for a G (0,2), 



/ Wr.R{dz) 




Je 





L2K 



(a-i)r°-i for a e (1,2) 
L2K for a e (0, 1). 



4. Convergence estimates for the discrete-time problem 

In this section we prove a priori error bounds for the convergence of solutions 
Uh of the semi-discrete scheme (I3.8P to the unique viscosity solution u of (|l.ip . We 
consider two cases: (i) The measure v is bounded and (Bl)o holds and (ii) the 
measure v is the truncation of an unbounded measure satisfying (B3). 

In view of the equi-boundedness, equi-continuity and consistency results of Propo- 
sitions EH] and [231 '^h converge locally uniformly to u by the Arzela-Ascoli Theorem, 
stability and uniqueness results for viscosity solutions (see e.g. [4j), and the con- 
sistency result in Proposition 13. II It is also possible to obtain convergence without 
equi-continuity (i.e. under weaker assumptions on the coefficients) using so-called 
half relaxed limits [7]. Such results are given in 14J for some non-local equations, 
but these results does not cover the HJB equations we consider here. 

Now we proceed to obtain a priori estimates for the convergence of Uh to u. To 
do this we will make use of an abstract result in ,,30] . which we will describe below. 
Consider the equation 

(4.1) F{x, u, Du, D^u, u(-)) = x G 

where u(-) represents non-local (integral) terms. Let ft, > be a discretization 
parameter and consider an approximation scheme for (j4.ip written in abstract form 
as 

(4.2) S{Kx,Uh{x),\uh\^)^Q xeM^, 

where [u^Ja, represents a function defined at x via all the possible value of u^. We 
need the following set of assumptions. 

(CI) (Monotonicity) There exists co > such that for any h > 0, x ^ R^, C G I^, 
r > and bounded functions u, v such that u < i> in R^, then 

S{h,x,( + T, [u + t]^) > S{h,xX, [v]x) + cqt 

(C2) (Regularity) For any h > and any continuous, bounded function the 
function 

X i-> S{h, X, 4){x), [4)]x) 
is bounded and continuous on R^ and the function 

C ^ S{h,x,C, [(f>]x) 
is uniformly continuous for bounded uniformly in x. 

(C3) ( Consistency) There exists a function E{K, h, e) such that for any sequence 
{0e}e of smooth functions satisfying 

\D^(t>,{x)\ < i^e^^l^l in R^, for any f3 e N^, 

where = X^iLi following inequality holds: 

\Sih, X, Ux), [0.].) - F{x, 0, i?2^,)| < E{k, h, e) in R^. 

(C4) (Convexity) Let (pe)£>o be a family of moUifiers (smooth, positive func- 
tions with mass 1 and support in {\x\ < e}). For any Lipschitz-continuous 
function 0, there exists a constant C such that for any x and h 

/ S{h,x,(t){x - e), [(/)(• - e)]x)pie)de > S{h,x, (0 * Pe)ix), [(</> * ps)]x) - Ce 
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(C5) (Commutation with translations) For any h > small enough, < e < 1, 
y e M^, C g M, continuous bounded function (j) and |e| < y, we have 

S{h, y, C, [(t)]y-e) = S{h, y, C, [(/>(• - e)]^^). 

(D) For h small enough and e E [0,1), there is a unique solution w| of the 
scheme 

(4.3) max5(/i,a; + e,u|(x), [u|]^) = in R^, 

\e\<e 

where Uh := w° also solve (|4.2p . and a constant C independent of h,e such 
that 

\ulU<C and |<-<|o<C£. 

We remark that we are using a more general consistency relation here than in 
[30j . and that this extra generality will be needed when we consider unbounded 
measures v. The next result is a restatement of Theorem 3.4 in [30] in view of the 
new consistency relation (C3). 

Theorem 4.1. Assume (Al) - (A3), (Bl)2, (CI) - (C5), and (D) hold, and let u 
and Ufi be solutions of respectively (j4.ip and (I4.2p satisfying K := \u\i V < oo. 
Then there exists a constant C depending only on Li,L2,co from (A2) and (A3) 
such that 

\u-Uh\<Cxmn{e + Ei{k,h,t)^ in R^. 

Remark 4.2. To prove (D) for the scheme p.Sp . we will need to assume also (Bl)o, 
h\ < Cq, h < 1, and cq > ci for both ci defined in Theorem 12.11 and Proposition 
13.31 Under these assumptions we also have K := V \uh\i < oo. 

We will apply this abstract result to derive error bounds for the scheme 
We rewrite the scheme in the form (14.21) with 



S{h,x,r, [i)]^) = 
(4.4) 

r -e"^'* r 

\ ojt. [■ilj]x{hb{x, v) + Vha„i{x, v)) - 2r + [4;]x{hb{x, v) ~ Vh(T,n{x,v)) 

^gy I. Zan — ' L 

m— 1 

1 _ p->^fi f phc{x.v) _ 1 ,^ 

^ + ry(x, V, z)) - r] + r - e^'^^^^^^ f{x, v)], 

and ['0]a;(2:) — Tpix + z). 

Lemma 4.3. Assume (Al) - (A3) and (Bl)o, hX < Cq, h < 1. Then the scheme 
p.8|) (and equivalently (|4.4p ) satisfies assumptions (C1)-(C5) with 

E{K,h,e)^Cih{Ke-^ +Ke-^ +Ke-^) + C2hX{{l + \ [ t]i^\)K + Ke-^), 

Je 

where the constant Ci and C2 only depend on sup„ |cr|o,sup„ |fe|o. 

// in addition cq > ci for both ci 's in Theorem \2.1\ and Proposition \ 3.!A then 
assumption (D) also holds with constants C only depending on the data and Cq. 

Proof. It is straightforward to verify (CI) with cq = cq where cq is defined in 
(A3). (C2) follows from the assumptions on the coefficients, while (C3) follows 
from Proposition 13.11 (v). By a straight forward computation, it follows that (C4) 
holds with C — 0. We refer to for similar computations. Finally, (C5) holds 
since [^Ji-e = [0(' — e)]x for any continuous function 0. 

To prove (D), observe that (|4.3p can be rewritten in the form (|4.4p (by defining 
a new control v — (w,e)). The coefficients of this new equation still satisfy (Al) - 
(A3). Therefore (D) follows after an application of Propositions 13.11 and 13.31 □ 
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Now we are in a position to state the error bounds. First we consider the bounded 
case, i.e (Bl)o holds. In this case the equation is (jl.ip which is approximated by 
the scheme (|3.8p . Note that the integral operator has not yet been discretized. 

Theorem 4.4 (Bounded measure). Assume (Al) - (A3), (B1)q, h<l, andcQ > c\ 
for both ci 's in Theorem \2.1\ and Provosition \3.3[ Let u be the solution of (jl.ip and 
Uh be the solution of (|3.8p (or equivalently ()4.4p ). 

(a) (General IPDEs) Then \u ~ Uh\o < Ch^/^. 

(b) (1st order IPDEs) If in addition a = 0, then \u — Uh\o < Ch^^^. 
In both cases the constant C depends only on the coefficients, ci, and A. 

Proof. Part (a) is an easy consequence of Theorem 14.11 and Lemma 14.31 Part (b) 
follows in a similar manner after noting that the consistency relation corresponding 
to Proposition [5Tl] (v) now becomes 

|LXx) - Ll^{x)\ < Cih\D^(t>\o + C2hX{{l + \ f w\)\D^\o + \D^^\o). 

J E 

□ 

Remark 4.5. The convergence rate obtained in Theorem 14.41 is the same as in the 
pure PDE case, see [27l [5]. In the first order case and when the Levy measure is 
bounded, convergence rate 1/2 for semi-Lagrangian schemes like (|3.8p have previ- 
ously be obtained in |16| . However the integral term in ^16) has a different form 
compared to the one we consider here. 

Remark 4.6. The scheme p.Sp uses a first order accurate approximation of 2nd 
derivatives as can be seen from the consistency relation Proposition 13.11 (v). This 
leads to lower rates of convergence than for some monotone FDMs that use 2nd 
order accurate approximations of 2nd derivatives, see [33l [30l [10] . There the rate is 
1/2, while for a more general class of FDMs the rate is at least 1/5, see for the 
pure PDE case. 

We proceed to the case of general unbounded Levy measure u under assumption 
(B3). There are two different cases: (i) a G [0,1) with HJB equation (jl.ip . and 
(ii) a e [1, 2) with HJB equation (|1.2p . In case (i), a G (0, 1), as a consequence of 
Lemma 331 Theorems 14.11 and 12.31 we have the following convergence result. 

Theorem 4.7 (Unbounded measure I). Assume (Al) - (A3), (B3) with a G (0, 1), 
h < 1, and cq > ci for both ci 's in Theorem \2.1\ and Provosition \3.3l Let u be the 
solution of (jl.ip and Uh be the solution of (I3.16p . 

3 

Then the best rate is obtained choosing r — /je+a , and in this case 

\u - Uh\o < C{h^^^ + h^) < Ch^l^, 

where the constant C depend only on the coefficients, Ci, and quantities from (B3)/(B1)2. 

Proof. When a G (0,1), Lemma [4.31 still holds for the scheme p.l6p . and in view 
of l|3.19p we have the following form of £', 

E{K, h, r, e) = CiKh{e-^ + e'^ + e"^) + C2Khr-"{l + e^^), 

where the constant Ci and C2 are independent of h, r, s. Let Ur denote the solution 
of (|2.2p . Theorem 14 . 1 1 and ( "term- wise" ) minimization in e, lead to the bound 

\ur - Uh\a < C{h^^^ + r-^h + r^^/i^), 

where the constant C depend only on the coefficients, ci, and quantities from 
(B3)/(B1)2. In view of Theorem 12.31 and the optimal choice of r, r — HbT^ , the 
result follows. □ 
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Remark 4.8. Since j < < i for a e (0, 1), there is no reduction of rate due to 
truncation of the measure v. 

The case a e (1, 2) is more difficult, since now also riVr,R and hence 6 in p.l7p 
blows up as r — > 0. As a consequence Theorem 14 . II can no longer be used directly. 
The convergence result is the following: 

Theorem 4.9 (Unbounded measure II) . Assume (Al) ~ (A3), (B3) witha e (1,2), 
h < I, and cq > c\ for both ci 's in Theorem \2.1\ and Proposition \3.3l Let u he the 
solution of (|1.2p and Uh he the solution of p.l7p . 

Then the best rate is obtained choosing r = h'^+^ , and in this case 

\u - Uh\o < Ch^, 

where the constant C depend only on the coefficients, Ci, and quantities from (B3)/(B1)2. 
Remark 4.10. This result is consistent with Theorem 14.71 since lim — h- For 

' ' -I I 3+5a 4 

a e (1, 2) the rate degrades as a increases, and lim i^^r" = t^-- 



Outline of proof. From Proposition 13.31 we have uniform in r Lipschitz continuity, 
but the continuous dependence estimates will be proportional to Jg \z\vr^R through 
the 6-term. Because of this we must redo the arguments leading to Theorem 14. 1[ 
and the result will be an estimate of the form 



Uh\<C min ( e / \z\vr,R. + Ei{K, h, r, e) 
^>o V Je 



when Ur solve (|2.6p . We omit the details, since the argument is same as the one 
used to prove Theorem 14. II in [30] . Because of the blow up in Wr,R, we also need 
a much more precise consistency relation than given in Proposition 13. II (v) tracking 
all b dependence. From the proof of Proposition 13.11 it is easy to see that it will 
have the following form 



E{K,h,e) = CiK h\b\le-^ + {h^\b\l + h\b\a)e-^ + ih + h^\b\^ + h^\b\l)e-'^ 



C2KXrh 



\b\o + e- 



+ C^KXrh 



Wr,R 

E 



where the constants Ci, C2, C3 only depend on sup„ |(t|o. In view of p.l9p . = 
Cr~" and |&|o+ l^^l^^r,/? < C(H-r"^^). The rest of the proof is a long computation 
consisting of choosing optimal e and r as in the proof of Theorem 14. 71 We omit the 
details only remarking that the worst term in E turns out to be the ft.e^'^-term. □ 

Appendix A. Lipschitz regularity and continuous dependence 

In this section we obtain a combined Lipschitz regularity and continuous depen- 
dence result for solutions to the equation 

(A.l) 



Uhix) = inf \hf{x, v) + e '"''-^''"^ V (uh{x + hb{x, v) + Vham{x, v))+ 

vev I. I za ^ — ' V 

m— 1 

+ Uh{x + hb{x,v) - Vham{x,v)^ + - — ^ / Uh{x + r]{x,v, z))i^{dz) | 



where A — v{E) and 



E 
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Theorem A.l. Let uy^ and he solutions of (jA.ip corresponding to the data 
cr, 5, c, /, 77, v and a, b, c, /, 77, v respectively, and assume both sets of coefficients sat- 
isfy (Al) - (A3) and (B1)q, and that hX < C and h < I. Then there exist constants 
ci,L,K > Q (only depending on the data and C) such that if cq > ci, then for all 
h>0,x,y& M^, 

\uh{x) - uh(y)\ < L\x - y\ +Ks\ip [\f - f\o + |c - c|o 

vev 

+ \b- b\o + W-a\o + \J H; z) - V{; z)\Mdz)\l^^] ■ 

Remark A. 2. The precise dependence of the constants ci, L, K is given in the proof 
below. 

Remark A. 3. This resuh extends the corresponding resuhs of [5] to non-local HJB 
equations, with general (singular) Levy measures. Moreover, the proof below sim- 
plifies the corresponding proofs of Barles and Jakobsen [5] because we do not use 
the somewhat unnatural "doubling schemes" as in [5] . Instead we give a more direct 
proof. 

Proof. We will use doubling of variables techniques similar to those used to prove 
corresponding results for equation p.ip . We define 

(A.2) y)^m + aM + L{a-^ + a\x - yf) + e{\xf + \yf), 

(A.3) ip{x, y) = sup[u,,(a;) - Uh{y) ~ (j){x, y)] = tp{x, y), 

x,y 

where a, £ > 0, m, M, L > 0, and (x, y) is the point where the supremum is attained. 
We will prove that ip{x,y) < o(l) as e ^ for a suitable constant L, and for 

(A.4) 

1/ - f\o + {\uh\o A \uh\o + h\f\o A |/|o)|c - c|o] , 
(A.5) 

M k[—^) sup [\a - a\l + \b~b\l + \j |r;(, z) - f,{; z)\^v{dz)\,] . 
This implies Theorem lA.il (after sending e ^ 0) because for any a e M and x,y E 

Uh{x) - Uh{y) - m - 3(LAf )^/^ - L\x - y\ < ip{x, y) < ip{x, y) < o(l) as e ^ 0. 

For simplicity we will not be explicit about the form of the e-terms appearing 
in the computations below. The role of these terms is only to guaranty that the 
maximum is attained at a point (x, y), and they vanish in the final estimate when 
e ^ 0. We refer to the proof of Theorem 3.4 in |5j for details concerning the e-terms. 

We proceed by contradiction assuming that ip{x,y) > o(l) as e ^ 0. Note that 
by the definition of ip this implies that Uh{x) — Uh{y) > 0. Furthermore, observe 
that since {x,y) is a maximum point, 

■0(2; + h + a,y + h + a) + ip{x + b ~ a,y + b — a) < 2ip{x, y), 

i(;{x + C,y + C) < ip{x,y), 

for every a, b, ^, a, 6, C G M.^ , and hence by the definition of tp, 

h '■= [uh{x + 6 + a) — 2vh{x) + Uh{x + 6 — a)] 

- [uh{y + b + a) - 2uh{y) + Uh{y + b-a)] 

< (j){x + b -\- a,y + b -\- a) — 2(j){x, y) -\- 4>{x + b — a,y + b — a), 

h ■= [uh{x + C) - Uh{x)] - [uh{y + C) - ^h{y)\ < <l>{x + C,y + C) - <l>{x, y)- 



sup e 



c|oV|c|o/i 
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Finally, by the definition of (j) we are lead to 

(A.6) Ii<2La\a-af + 2La\b-b\'^ +4La{x-y,b~b)+o{l) as e ^ 0, 
(A.7) h < 2La{x -y,C-C) + La\C - Cl' + 

These two inequalities are crucial for the rest of the proof. 
Now we divide (jA.ip by he~^'^ and rewrite it as 



as e ^ 0. 



= sup <^ Uh{x) - e ''f{x, v) 

vev ^ n 

2dh ^ ("'»(^ + hb{x,v) + Vha„i{x, v)) - 2uh{x) + Uh{x -\- hb{x,v) - \/ham[x, v) 



m—l 

1 - e-^^ 



(^Uh{x + ri{x,v, z)) - Uhix)ji^{dz)Y 



Upon subtracting the equations (in this new form) for Uh and Uh, we get 



< sup I 



-Uhiy) 



Uh{x) 



+ [e^>^f{x,v)-e''^f{y,v)] + -—J2h^ 



1-e 



2dh ^ ■ \h 

rn—1 



2v{dz)^, 



where Ii.m corresponds to Ii above with the choice a = \/ham{x,v), b = hb(x,v), 
a = y/ham{y,v), and b = hb{y,v), and for I2 we have taken C = r](x,v,z) and 
( = r]{y, V, z). By (jA.6P and (|A.7P we then have 

h 1 



< sup 



I — — Uh{y) — Uh{x) 



+ [e=V(S, v) - e^f^fiy, v)] + -^La [- ^ 2h\a„,{x, v) ~ v)\ 



+ 2h^\b{x, v) - b{y, + Ah{x - y, b{x, v) - b{y, v)) 



+ e 



-he 



-hX 



Xh 

+ 0(1) as £ — > 



La 



\r]{x,v, z) - fi{y,v,z)\^ + 2{x - y,-q{x,v,z) - ri{y,v, z)) iy{dz)^ 



Since b — b + \^-xh Je ^ind b is defined similarly we see that the {x — y,T] — f])- 
terms cancel in the above inequality. Since v{E) — A, Jensen's inequality implies 
that 



v{dz) 



{ri{x,v,z) -fj{y,v,z)) 



Also note that since Uh{x) — Uh{y) > and c > cq > 0, 



< AM \7]{x,v,z) -f}{y,v,z) 



, v{dz) 



-Uh{y)- 



Uhix) < 



'uhix)-Uhiy)]+ 



ch ^ch\ 



\uh{x)\A\uh{y)\. 
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Therefore after cancellations, Jensen's inequality, and the inequality < he'^' 

we get 



[uh{x) - Uh{y)] 



< sup 



|e|c|ov|c|o^ v) - f{y, v)\ + {\u^\o A \uh\o + h\f\o A |/|o)|c(x, v) - c{y, v)\ 

-h\ ^ d 



m— 1 

+ Ah^\b{x, v) ~ b{y, w)p + Ah{x - y, h{x, v) - b{y, v)) 



1 - e-^'^ ^ 



Ah'-'Xhe^'' I \r]{x,v,z) -fi{y,v,z)\'^i^{dz) 



+ \ri{x,v, z) — fi{y,v, z)\'^ i^{dz)^ + o{l) as e — > 0. 
Now since Xh < C and h < 1, it follows from simple computations that 

7 ['0(^,2/) +Hx,y)] = [uh{x) - Uhiy)] 



< sup jel^l"''!'^!"'' 



h 

1/ - /|o + i\uh\o A \uh\o + h\f\o A |/|o)|c - c|o 



m—1 



LahCe^\ 



\'q{-,v,z) -r)(-,i;,z)pjy(dz)|o| 
+ sup jel^l"^!'^!"'' \Lf + {\uh\Q A |u,,|o + /i|/|o A |/|o)icl + a\x - y\) 



La\x-y\\^YdY.^^l + ^^hL 



6 + 8 



+ o(l) as e ^ 0. 
Let L, Lq be defined as in the inequality above. If cq is so big that 



^cah _ 2 



Lo > 0, 



and we choose m, M as in (|A.4p and (jA.Sp for K big enough, and 



-in 



then V'la;, y) < o(l) as e ^ and the proof is complete. 

References 



□ 



[1] A. Almendral and C. W. Oosterlee Accurate evaluation of European and American options 
under the CGMY process. SIAM J. Sci. Comput. 29 (2007), no. 1, 93-117 

[2] S. Asmussen and J Rosinski. Approximations of small jumps of Levy processes with a view 
towards simulation. J. Appl. Probab. 38(2): 482-493, 2001. 

[3] G. Barles, R. Buckdahn, and E. Pardoux. Backward stochastic differential equations and 
integral-partial differential equations. Stochastics Stochastics Rep., 60(l-2):57-83, 1997. 

[4] G. Barles and C. Imbert. Second-Order Elliptic Integro-Differential Equations: Viscosity So- 
lutions' Theory Revisited. Ann. Inst. H. Poincare Anal. Non Lineaire, 2007, to appear. 



FINITE ELEMENT SCHEME FOR INTEGRO-PARTIAL DIFFERENTIAL EQUATIONS 21 



G. Baxles and E. R. JaJjobsen. On the convergence rate of approximation schemes for 
Hamilton- Jacobi-Bellman equations. M2AN Math. Model. Numer. Anal. 36(1): 33-54, 2002. 
G. Barles and E. R. Jakobsen. Error bounds for monotone approximation schemes for 
Hamilton-Jacobi-Bellman equations. SIAM J. Numer. Anal. 43(2):540-558, 2005. 
G. Barles and P. E. Souganidis. Convergence of approximation schemes for fully nonlinear 
second order equations. Asymptotic Anal. 4(3): 271-283, 1991. 

F. E. Bentli, K. H. Karlsen, and K. Rcikvam. Portfolio optimization in a Levy market with 
intertemporal substitution and transaction costs. Stock. Stock. Bcp., 74(3-4):517-569, 2002. 

D. P.Bertsekas and S.Shreve. Stochastic optimal control: the discrete time case. Academic 
Press, New York, 1978. 

I. H. Biswas, E. R. Jakobsen, and K. H. Karlsen. Error estimates for finite difference- 
quadrature schemes for fully nonlinear degenerate parabolic integro-PDEs. J. Hyperbolic Dif- 
fer. Equ., 5(1): 187-219, 2008. 

I. H. Biswas, E. R. Jakobsen, and K. H. Karlsen. Viscosity solutions for a system of integro- 
PDEs and connections to optimal switching and control of jump-diffusion processes. Submit- 
ted, 2008. 

F. Bonnans and H. Zidani. Consistency of generalized finite difference schemes for the sto- 
chastic HJB equation. SIAM J. Numer. Anal. 41(3);1008-1021, 2003. 

S. C. Brenner and L. R. Scott. The mathematical theory of finite element methods. Springer- 
Verlag, New York, 2002. 

M.Briani, C. La Chioma, and R.Natalini. Convergence of numerical schemes for viscosity solu- 
tions to integro-differential degenerate parabolic problems arising in financial theory. Numer. 
Math. 98 (4):607-646, 2004. 

M. Briani, R. Natalini, and G. Russo. Implicit-explicit numerical schemes for jump-diffusion 

processes. CALCOLO, 44(1): 33-57, 2007. 

F. Camilli Approximation of integro-differential equations associated with piecewise deter- 
ministic process. Optimal Control Appl. Methods 18(6): 423-444, 1997. 

F. Camilli and M. Falcone, An approximation scheme for the optimal control of diffusion 
processes, RAIRO Model. Math. Anal. Numer. 29:97-122, 1995. 

P. Carr, H. Geman, D. B. Madan, and M. Yor. Stochastic volatility for Levy processes. Math. 
Finance 13 (3):345-382, 2003. 

R. Cont and P. Tankov. Financial modelling witk jump processes. Chapman & Hall/CRC 
Financial Mathematics Series. Chapman & Hall/CRC, Boca Raton, FL, 2004. 
R. Cont and E. Voltchkova. Finite difference methods for option pricing in jump-diffusion and 
exponential Levy models. SIAM J. Numer. Anal. 43(4): 1596-1626, 2005. 
M. G. Crandall and P.-L. Lions Two approximations of solutions of Hamilton-Jacobi equa- 
tions. Math. Comp. 43(167): 1-19, 1984. 

S. Elganjoui and K. H. Karlsen. A markov chain approximation scheme for a singu- 
lar investment-consumption problem with Levy driven stock prices, available at the URL 
http : //www . math . uio . no/ "kennethk/ , 2002. 

Y. d'Halluin, P. A. Forsyth, and G. Labahn. A penalty method for American options with 

jump-diffusion processes. Numer. Matk., 97(2):321-352, 2004. 

M. Falcone Appendix A in M. Bardi and I. Capuzzo Dolcetta Optimal Control and Viscosity 
Solutions of Hamilton-Jacobi-Bellman equations, Birkhauser, Boston, 1997. 
W. H. Fleming and H. M. Soner. Controlled Markov processes and viscosity solutions. Second 
edition. Springer, New York, 2006. 

I. I. Gihman and A. V. Skorohod Stochastic differential equations. Translated from the Rus- 
sian by Kenneth Wickwire. Ergebnisse der Mathematik und ihrer Grenzgebiete, Band 72. 
Springer- Verlag, New York-Heidelberg, 1972. 

E. R. Jakobsen. On the rate of convergence of approximation schemes for Bellman equations 
associated with optimal stopping time problems. Math. Models Methods Appl. Sci. (M3AS) 

13(5): 613-644, 2003. 

E. R. Jakobsen and K. H. Karlsen. Continuous dependence estimates for viscosity solutions 
of integro-PDEs. J. Differential Equations, 212(2):278-318, 2005. 

E. R. Jakobsen and K. H. Karlsen. A "maximum principle for semicontirmous functions" 
applicable to integro-partial differential equations. NoDEA Nonlinear Differential Equations 
Appl. 13:137-165, 2006. 

E. R. JaJfobsen, K. H. Karlsen, and C. La Chioma, Error estimates for approximate solutions 
to Bellman equations associated with controlled jump-diffusions. Submitted 2005. 
N. V. Krylov. On the rate of convergence of finite-difference approximations for Bellman's 
equations. Algebra i Analiz, 9(3):245-256, 1997. 

N. V. Krylov. On the rate of convergence of finite-difference approximations for Bellman's 
equations with variable coefficients. Probab. Theory Related Fields, 117(1):1-16, 2000. 



22 



FABIO CAMILLI AND ESPEN R. JAKOBSEN 



[33] N. V. Krylov. The rate of convergence of finite-difference approximations for Bellman equa- 
tions with Lipschitz coefficients. Appl. Math. Optim. 52(3):365-399, 2005. 

[34] H. J. Kushner and P. Dupuis. Numerical methods for stochastic control problems in continuous 
time. Second edition. Applications of Mathematics 24. Stochastic Modelling and Applied 
Probability. Springer- Verlag, New York, 2001. 

[35] P.-L. Lions Existence results for first-order Hamilton-Jacobi equations. Ricerche Mat. 32(1):3- 
23, 1983. 

[36] A.-M. Matache, T. von PctcrsdorfT, and C. Schwab. Fast deterministic pricing of options on 
Levy driven assets. M2AN Math. Model. Numer. Anal, 38(1):37-71, 2004. 

[37] J. L. Menaldi Some estimates for finite difference approximations. SIAM J. Control Optim. 
27: 579-607, 1989. 

[38] R. Merton. Option prices when the underlying stock returns are discontinuous. J. Financial 

Economics, 3:125 144, 1976. 

[39] H. Pliaiu. Optimal stopping of controlled jump- diffusion processes: a viscosity solution ap- 
proach. J. Math. Systems Estim. Control, 8(1), 1998. 

[40] T. H. Rydberg. The normal inverse Gaussian Levy process: simulation and approximation. 
Heavy tails and highly volatile phenomena. Comm. Statist. Stochastic Models 13(4):887- 
910,1997. 

[41] A. Sayah. Equations d'Hamilton-Jacobi du premier ordre avec termes integro-difierentiels. L 
Unicite des solutions de viscosite. Comm. Partial Differential Equations, 16(6-7):1057-1074, 
1991. 

Dip. di Matematica Pura e Applicata, Univ. dell'Aquila, log. Monteluco di Roio, 
67040 l'Aquila, ITALY 

E-mail address: camilli®ing. univaq.it 

Department of Mathematical Sciences, Norwegian University of Science and Tech- 
nology. 7491 Trondheim, NORWAY 
E-mail address: erj@matli.ntnu.no 



